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Qh' The muon charge ratio of ultrahigh energy cosmic rays may provide information to 

detect the composition of the primary cosmic rays. We propose to extract the charge 
information of high energy muons in very inclined extensive air showers by analyzing 
their relative lateral positions in the shower transverse plane. 



CZ3 

The most high energy particles can be observed by human being are from cosmic 
rays. The study of them belongs to frontiers of human knowledge in combination of 
, cosmology, astrophysics, and particle physics, and can provide better understanding 

qq ' of the universe from most small to most big, i.e., connecting quarks to the cosmos. 

■ The universe is not empty, but full of background relic particles from the big bang. 

It has long been anticipated that the highest energy cosmic rays would be protons 
from outside the galaxy, and there is an upper limit of the highest energy in the 
observed proton spectrum, commonly referred to as the GZK cutoff - , as the protons 
| traveling from intergalactic distances should experience energy losses owing to pion 

productions by the photons in the cosmic background radiation. Although there 
have been attentions for the cosmic ray events above the GZK cutoff, it is natural 
to expect that these ultrahigh energy cosmic rays come from sources within the GZK 
zone 2, i.e., not far from us in more than tens of Mpc. Recentlythere are also reports 

To] 

on the observation of the GZK cut-off by new experiments Q . However, questions 
about the composition of such ultrahigh energy cosmic ray particles, e.g., whether 
they are protons, neutrons, or anti-nucleons are still open to investigations. 

Muons in the air showers are mainly from decays of pions and kaons produced 
in the interactions of the primary cosmic rays with the atmosphere. The very high 
energy secondary pion and kaon cosmic rays can be considered as from the current 
fragmentation of partons in deep inelastic scattering of the primary cosmic rays with 
the nucleon targets of the atmosphere in a first approximation^. We also consider 
only the favored fragmentation processes, i.e., the 7r + , which is composed of valence 
u and d quarks, is from the fragmentation of u and d quarks in the nucleon beam, 
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and the ir~ , which is composed of valence u and d quarks, is from the fragmentation 
of u and d quarks^!. Similarly, the K + , which is composed of valence u and s, is from 
the fragmentation of u and s quarks, and the K~ , which is composed of valence 
u and s, is from the fragmentation of u and s quarks. The fi + is from the decay 
of a 7r + or a K + and the /i~ is from the decay of a tt~ or a K~ . We can roughly 
estimate the muon charge ratio by 

/i+ = Jq dx { [u(x) + d(x)] + ac [u(g) + s(x)} } 

M~ da; {[d{x) + u(a;)] + as [u{x) + s(x)]} ' 

where q(x) is the quark distribution with flavor q for the incident hadron beam 
and k ~ 0.1 — » 0.3 is a factor reflecting the relative muon flux and fragmentation 
behavior of K/tt. Secondary collisions do not influence the above estimation, since 
the current parton beams still keep their flavor content and act as the current 
partons after the strong interactions with the partons in the atmosphere targets. 
Adopting a simple model estimation of the parton flavor content in the nucleon 
without any parameter^, we find that ~ 1.7 for proton and ~ 0.7 

for neutron. This simple evaluation is in agreement with the empirical expectation 
of 1 1 [sT k, 1.66 for proton and fi + / fjT w 0.695 for neutron^ as well as that in an 
extensive Monte Carlo calculation^, thus it provides a clear picture to understand 
the dominant features for the muon charge ratio by the primary hadronic cosmic 
rays. For the ratio for antiproton, it is equivalent to the /z - //i + ratio for 

proton by using Eq. (fTJ), thus we find [i + / fi~ ~ 0.6 for antiproton, which is close to 
that for neutron. The fi + / /i~ ratio for antineutron is also equivalent to the 
ratio for neutron, and it is ~ 1.4, which is close to that for proton. It 

is hard to distinguish between the primary neutrons and antiprotons (or protons 
and antineutrons) by the ratio of the air shower, unless very high precision 

measurement is performed and also our knowledge of the muon charge ratio for each 
nucleon species is well established. 

The study of cosmic rays with primary energies above 10 5 GeV are typically 
based on the measurements of extensive air showers (EAS) that they initiate in the 
atmosphere. The ground detector array records the secondary particles produced in 
shower cascades, including photons, electrons (positrons), muons, and some hadrons. 
Then their arrival times and density profiles are used to infer the primary energy 
and composition of the incident cosmic ray particle, usually through comparison 
with simulated results. Photons, electrons and positrons are the most numerous 
secondary particles in an EAS event. However, for very inclined showers, these 
electromagnetic components would travel a long slant distance and are almost com- 
pletely absorbed before they reach the ground. On the other hand, muons are decay 
products of charged mesons in shower hadronic cascades. Most high energy muons 
survive their propagation through the slant atmospheric depth, during which they 
lose typically a few tens of GeV's energy. These high energy muons carry impor- 
tant information about the nature of the primary cosmic ray hadron, which will be 
extracted from their energy spectrum and lateral distribution. 
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As discussed by Hwang and I ^, the ratio of positive versus negative muons 
H + / fi~ is a significant quantity which can help to discern the primary composition, 
and at high energies this charge ratio also reflects important features of hadronic 
meson production in cosmic ray collisions. In order to obtain such muon charge 
information, we would need a way to distinguish between positive and negative high 
energy muons. Unfortunately, existing muon detectors available at shower arrays, 
usually scintillators and water Cerenkov detectors, are not commonly equipped with 
magnetized steel to differentiate the muon charges. Even if they were, the limited 
region of the magnetic field prevents definite determination of high energy muons' 
track curvature. 

This invites us to think of the geomagnetic field as a huge natural detector for 
muon charge information. Apparently, after being produced high in the atmosphere, 
a positively charged muon would bend east on its way down while a negatively 
charged muon would bend west, introducing an asymmetry into the density profile 
of the shower front. If their separation is large enough as compared with other 
circularly symmetric "background" deviations, it will be possible to distinguish the 
positive muons from the negative one s. 

To see such an effect, Xue and I analyzed the possibility of obtaining the 
charge information of high energy muons in very inclined extensive air showers. 
We have demonstrated that positive and negative high energy muons in sufficiently 
inclined air showers can be distinguished from each other through their opposite ge- 
omagnetic deviations in the transverse plane. We developed a revised Heitler model 
to calculate this distinct double-lobed distribution, and studied the condition for 
the two lobes of either positive or negative muons to be separable with confidence. 
From our criterion of resolvability, we concluded that a zenith angle 75° < 9 < 85° 
will be most suitable for our approach. 

There are already some results from full air sho wer simulations that tak e into 
account the geomagnetic effect on muon propagation tllll211311411all6J17J18J xhey 
illustrated remarkable double-lobed muon lateral density profile in very inclined 
air showers, which is in agreement with our expectation qualitatively. However, no 
present study has fully considered the high energy part of muon content, which 
can be used to compare with our results. Thus we would like to propose future 
simulations of very inclined extensive air showers that focus on the behavior of high 
energy muons. They also have to keep track of the muon charges and the relation 
to their lateral positions. For more detailed analysis and discussion, please refer to 
RefP. 

In summary, we propose to extract the charge information of high energy muons 
in very inclined extensive air showers by analyzing their relative lateral positions in 
the shower transverse plane. This muon charge information is helpful to detect the 
composition of cosmic rays, e.g., the neutron or antiproton content of the ultrahigh 
energy cosmic rays. 



August 19, 2008 18:16 WSPC/INSTRUCTION FILE mabq 



4 BO-QIANG MA 
Acknowledgments 

I am very grateful to Pauchy Hwang and the organizers for their invitation and 
warm hospitality. I also thank Pauchy Hwang and BingKan Xue for the collab- 
orated results in this talk. This work is partially supported by National Natural 
Science Foundation of China (Nos. 10721063, 10575003, 10528510), by the Key 
Grant Project of Chinese Ministry of Education (No. 305001), and by the Research 
Fund for the Doctoral Program of Higher Education (China). 

References 

1. K. Greisen, Phys. Rev. Lett. 16, 748 (1966); G.T. Zatsepin and V.A. Kuzmin, Pis'ma 
Zh. Eksp. Tero. Fiz. 4, 114 (1966) [JETP Lett. 4, 78 (1966)]. 

2. F.W. Stecker, Phys. Rev. Lett. 21, 1016 (1968). For an extensive review, see, 
F.W. Stecker, astro-ph/0101072| 

3. R. Abbasi et al, HiRes Collaboration, Phys. Rev. Lett. 100, 101101 (2008). Also Pierre 
Auger Observatory reported the observation of the GZK cut-off. 

4. W-Y.P. Hwang and B.-Q. Ma, Eur. Phys. J. A 25, 467 (2005). 

5. See, e.g., B.-Q. Ma, I. Schmidt, J. Soffer, and J.-J. Yang, Nucl. Phys. A 703, 346 (2002). 

6. See, e.g., B.-Q. Ma, I. Schmidt, and J.-J. Yang, Phys. Rev. D 65, 034010 (2002). 

7. Y.-J. Zhang, B. Zhang, and B.-Q. Ma, Phys. Lett. B 523, 260 (2001). 

8. R.K. Adair, Phys. Rev. Lett. 33, 115 (1974); R.K. Adair et al, Phys. Rev. Lett. 39, 
112 (1977); O.C. Allkofer et al, Phys. Rev. Lett. 41, 832 (1978). 

9. J.N. Capdevielle and Y. Muraki, Astropart. Phys. 11, 335 (1999). See, e.g., Fig. 7. 

10. B. Xue and B.-Q. Ma, Astropart. Phys. 27, 286 (2007). 

11. M. Ave, R. A. Vazquez and E. Zas, Astropart. Phys. 14, 91 (2000). 

12. H. H. Aly, M. F. Kaplon, and M. L. Shen, Nuovo. Cim. 31, 905 (1964). 

13. J. R. Horandel, J. Phys. G 29, 2439 (2003). 

14. C. A. Ayre et al, J. Phys. A 5, L102 (1972). 

15. T. Hebbeker and C. Timmermans, Astropart. Phys. 18, 107 (2002). 

16. S. Ostapchenko, Nucl. Phys. Proc. Suppl. 151, 147 (2006). 

17. See, e.g., W. D. Apel et al. [KASCADE Collaboration] , |arXiv:astro-ph/05108io| 

18. P. Hansen, T. K. Gaisser, T. Stanev and S. J. Sciutt o, Phys. Rev. D 71, 083012 (2005); 
A. Cillis and S. J. Sciutto, |arXiv:astro-ph/9908002| 



